Koshino- Taylor effect in graphene 
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We discuss the phonon-assisted scattering of electrons by defects, i.e., the so-called Koshino- 
Taylor effect, in graphene. The two-dimensional character of graphene implies that the strength of 
the Koshino- Taylor effect can be considerably larger than in ordinary metals. We show that at finite 
temperatures the defect-induced resistivity formally diverges in the thermodynamic limit, having a 
non-analytic TlnT component when finite size effects are taken into account. 



The possibility of fabricating free-standing graphene 
sheets has recently been demonstrated experimentally 
[l|, [2| . In this form graphene combines its unusual elec- 
tronic features with additional mechanical peculiarities, 
giving rise to new specific properties (see, e.g., Ref. 
[3|). The resistivity of suspended graphene, for exam- 
ple, shows a clear dependence on temperature 0, 0| in 
contrast to samples supported by a substrate where it 
has an important T- independent contribution Q. 

The most distinctive feature of graphene is the Dirac 
spectrum of its charge carriers Q. Suspended graphene, 
in its own turn, is inevitably subjected to both in-plane 
and out-of-planc distortions. As shown in Ref. 
this gives rise to two different contributions in its low- 
temperature resistivity. Whereas the scattering of elec- 
trons with in-plane phonons provides a contribution 
(see also Q), the scattering with out-of-plane ones gives 
an unusual T^/^lnT behavior. On the other hand, the 
precise nature of the defects responsible for the resid- 
ual (temperature-independent) resistivity is still unclear 
in graphene 0, @1 ■ If they could be described in terms 
of Dirac-delta- function potentials, the resulting resistiv- 
ity would not depend on the carrier density n. In the 
experiments, however, it turns out to be inversely pro- 
portional to n 0. This non-trivial behavior can be ex- 
plained in different ways: due to unscreened Coulomb 
potentials [l3, 111, frozen corrugations of the graphene 
sheet [lH , the scattering involving the midgap states that 
may be created by disorder [Tst . etc. The renormalization 
of the impurity scattering due to e.g. electron-electron 
interaction is another source of temperature-dependent 
contributions to the resistivity. Specifically in graphene, 
Friedel oscillations in the exchange field yield a linear- 
in-T behavior (which is absent in the case of Coulomb 
scatterers) In this paper we consider the phonon- 

assisted scattering of electrons by defects, i.e., the so- 
called Koshino-Taylor effect [H, M, [B, E, El. This 
effect is operative for any kind of defect and, in the case 
of graphene, is expected to be unusually large due to the 
two-dimensional character of the system. As we show 
below, both in-plane and out-of-plane phonons can yield 
a TlnT behavior of the low-temperature resistivity in 
graphene through the Koshino-Taylor effect. 



The key point in Koshino-Taylor effect are the local 
fluctuations of the system as experienced by defects. In 
fact, in real systems, the position of the defects is not 
fixed but fluctuates, what eventually translates into a 
temperature dependence of the defect-induced resistiv- 
ity di, [3, H 111 



19j. In ordinary metals this a rather 
small effect that, nevertheless, can be verified experimen- 
tally after some effort [2^. In the case of free-standing 
graphene, conversely, an unusually large effect can be 
anticipated because of the following. Within a first ap- 
proximation, it can be assumed that the changes in the 
positions of the defects coincide with the lattice fluctu- 
ations at the corresponding points. Thus, the residual 
resistivity obtained under the assumption of static impu- 
rities is effectively dressed by a Debye- Waller-like factor 
due to these fluctuations. Local fluctuations in a two- 
dimensional (2D) lattice formally diverge and, in accor- 
dance with the Mermin- Wagner theorem, such a diver- 
gence should destroy the crystalline order. In practice, 
however, long-wavclcngth fluctuations arc suppressed by 
finite size effects for example. This stabilizes the lat- 
tice, but local fluctuations may remain quite large. The 
Koshino-Taylor effect, benefiting from these large local 
fiuctuations, is thus expected to be amplified in free- 
standing graphene sheets. 

The simplest way to confirm this expectation is by 
means of the Boltzmann transport theory. This theory 
is appHcable in the case of gated and/or doped graphene 
with a finite (and tunable) carrier density (i.e. the so- 
called extrinsic graphene) [1, 21 1. Despite interference 
contributions are not captured within this approach, it 
suffices to reproduce correctly the qualitative behavior in 
our case [l8|. The resistivity then can be written as 
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(1) 



where vp is the Fermi velocity, v is the density of states at 
the Fermi level, and r is the transport relaxation time at 
the Fermi energy averaged over the angles (see, e.g., Ref. 
[l^l). In graphene = 2kp/{TrhvF), where the Fermi 
wavevector kp can be expressed in terms of the carrier 
concentration a,s hp = (tttt,)^/^. 

The relaxation time r is obtained from the so-called 
collision integral [l^l- So let us reconsider, first of all, 
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the scattering of an electron due to a single defect orig- 
inally situated at r = 0. The scattering potential asso- 
ciated with the defect is V{r — u), where u represents 
its displacement. Let k and k' be the initial and final 
wavevectors of the electron, and |A^) and |A^') the initial 
and final phonon states respectively. The scattering rate 
then can be written as 

-\V^',^\'\{N'\e^^--\N)\H{Ef - E,). (2) 

Here Vk'.k is the matrix element one would obtain by 
neglecting the displacement of the defect (i.e. taking u = 
0); K = k' — k is the change in the wavevector of the 
electron, and Ei {Ef) represents the total energy of the 
initial (final) state. It is worth noticing that the vector 
K lies on the graphene plane, so only the scattering with 
the in-plane modes turns out effective here. 

In the case of small displacements, it suffices to retain 
the lowest order terms of expansion of the exponential 



in powers of u: e'^ " ~ 1 + zK • u — ^(K • u)^. Fur- 
thermore, by assuming that the displacement of the de- 
fect coincides with the displacement of the lattice, i.e., 
the defect is pinned to the lattice, the vector u can 
be written as a sum of creation al" „ and annihilation 
operators Oq ^ for the phonon modes. Thus the main 
contribution to the scattering rate ([2]) is given by two 
types of processes. On one hand we have one-phonon 
processes, for which we have to consider the matrix ele- 
ments {N - l\a\N) = {N\a+\N - 1) = A^^/^ where N is 
the phonon distribution function. On the other hand we 
have processes involving no change in the phonon state, 
which are associated the matrix elements {N\N) = 1 
and {N\aa+\N) = {N\a+a\N) -1= N. Within our ap- 
proximations the graphene sheet can be considered as a 
continuum elastic medium with the two in-plane modes 
having the same velocity. Thus, after summation over 
the phonon states, the scattering rate for the electron is 
given by 
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r r hK^ n hK^ ~i 



(3) 



where a is the mass density and cjq = cq is the phonon frequency (here and hereafter the area of the system is taken 
as unity). 

In the case of a system of electrons, this scattering rate has to be multiplied by the electron distribution function 
/k for the initial state and by 1 — /k' for the final state. The collision integral then takes the form / = /o + /e/ + lim 
where 
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Xl^k',k|'(/k'-/k)<5(£k'-ek), 



hi = ~^n, J2 l^4',kp|^(/k' - /k)(l + 2iVq)5(£k' - £k) 



k',q 



2(TUJc, 



2tt hK^ 

hn = —n, |Vk'.kP^^{[./k'(l - /k)(l + ^q) + /k(l - /k')^q]'5(ek' - £k - nw^) 



k',q 



[/k'(l - /k)iVq + /k(l - /k')(l + A^q)]'5(£k' - £k + fi^q)}, 



(4a) 
(4b) 

(4c) 



with Ui being the density of defects. For the sake of 
concreteness, we assume in the following that the de- 
viation from equilibrium is due to the presence of a 
static and spatially homogeneous electric field. If such 
a deviation is small, Eqs. ([4]) can be linearized with 
respect to 5 = / — /□ = ~ /oIi^/q) ^^^ where 

/0(e) = (e^^^'')/^ + 1)^"'^, assuming that the phonon dis- 
tribution function takes its value of equilibrium: N = 

The integral (|la|) is nothing but the integral that one 



obtains neglecting the fluctuations in the positions of the 
defects. This further gives a residual resistivity pp, as 
can be seen in standard textbooks (e.g. Ref. [23|). As 
regards (j4bp , it can be written as lei = —0/ Tei after some 
standard manipulations (see e.g. Ref. [24I). Here 




where W{9) = ^ni|Vk',kP(l -cos6')2 {0 being the angle 
between k' and k) . At T = this gives a correction to the 
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residual resistivity obtained from the previous integral 
(|4a|) . This correction can be estimated as 



' MB 



Po, 



(6) 



where M is the mass of the unit cell and 8 is the Debye 
temperature. In the usual case ep = fi^kp/i^m), where 
m is the mass of the electron. In the above expression 
we then have -jj^ = jj^- This factor is typically ~ 
10~^ for ordinary metals. In the case of graphene ep ^ 
hvpkp, so the above factor can be written as 
This factor can reach the same value than in ordinary 



metals for large enough concentrations of charge carriers 
(i.e., large values of the Fermi energy ep). Actually at 
such a concentrations deviations from the Dirac spectrum 
might be important. However, as we can see, the precise 
form of the spectrum is not the crucial point for having 
the above correction to the residual resistivity due to the 
Koshino- Taylor effect. 



Let us now turn to the dependence on the temperature 
of the defect-induced resistivity. This dependence is due 
to both, the elastic processes giving rise to ([S]), and the 
inelastic processes taken into account in 



k'.q ^ 



[/o(ek') - /o(ek)]iV(fi£Oq)[l + N{hujq)]{^i,> - </?k)['5(ek' - £k - hw^) - (5(ek' 



(7) 
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In fact, the remaining contribution at low temperatures is 
due to this later integral (see Ref. [l3|)- This is basically 
because here we have the factors N{1 + N) whereas in 
(j4b|) we have N only. As a result of these factors the 
differences in the electron energies Ek' ^ £k = ztftwq are 
small, so we can put 



/o(£k') - /o(£k) ~ ±hhJq 



dh{e) 



(8) 



The integral ([7]) then takes the form w —g/Tim where 

Tin T-r^ crj 



k',q 



(9) 



The corresponding contribution to the resistivity can be 
estimated as 



h^klT (T^ 
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where Tg is temperature associated with the infrared cut- 
off of the integral over the phonon wavevectors. In the 
expression obtained for an ordinary (3D) metal one has 
an additional factor T/Q instead of ln(r/ro) O, [H]. 
This latter logarithmic factor is the fingerprint of the 
local fluctuations of the 2D lattice, whose largeness is 
ultimately determined by the quantity Tq (see below). 

Let us now discuss the possibility of a similar TlnT 
contribution due to the out-of-plane distortions of the 
graphene sheet. The virtual changes in the form of the 
defect potentials owing to the (local) deformations of the 
lattice brings about this possibility. That is, changes 
described as (1 -I- Xua)V{Y), where un is the trace of 



the strain tensor and A is a constant. This effect was 
first discussed in Ref. [2^ for ordinary metals. In 
the case of graphene the strain tensor takes the from 



djUi 



{dih){djhy\, where h character- 



izes the out-of-plane distortions [7[. Here we see that 
local in-plane deformations (uu) are possible as a result 
of out-of-plane fluctuations of the lattice. The resistivity 
then can be written as p ~ (1 -I- X{uii))po. In the absence 
of tension in the graphene sheet the dispersion law of the 
out-of-plane phonons is lu^ '^ = Kq^. At low temperatures 
we then have 



(uu) 



2ttMk 



T 



In 



T 



[h) 



(11) 



Here 0^'') is the Debye temperature for the out-of-plane 
phonons and T^^^ is associated with the infrared cut- 
off for these phonons (see below) . It is worth mentioning 
that if the graphene sheet is under tension, the dispersion 
law of the out-of-plane phonons becomes linear-in-g 0]. 
In this case the temperature dependence of {uu) changes 
from TlnT to , and the corresponding contribution 
to the resistivity is then similar to that obtained from 
the mean value of the square of the local deformations 
{ufi) (which is due to in-plane distortions). This sensitiv- 



ity to applied tension is absent in the contribution (fTOj) 
due to in-plane fluctuations. It is worth noticing also 
that these two contributions differ in their dependence 
on the Fermi wavevector kp, i.e., on the carrier density. 
In the out-of-plane contribution this dependence comes 
from the matrix elements of the (unchanged) scattering 
potential V only. Consequently, its dependence on the 
carrier density turns out to be the same than the one of 
bare resistivity po- This is not the case of contribution 
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due to in-plane phonons, where there is an additional 
dependence coming from the fluctuations [see Eq. pO|) ]. 

The dependence on the temperature of the resistivity 
is quite apparent in suspended graphene 0, Q ■ Far from 
the charge-neutrality point, this can be explained as due 
to electron-phonon scattering for temperatures above 50 
K 0, At lower temperatures this scattering is ex- 
pected to give and T^/^ In T contributions due to in- 
plane and out-of-plane phonons respectively So, if 
the temperature is low enough, the Koshino- Taylor ef- 
fect can overwhelm the electron-phonon scattering giv- 
ing rise to a T In T behavior in the resistivity. The range 
of temperatures at which the Koshino- Taylor effect can 
be dominant is difficult to estimate realistically given its 
strong dependence on the nature of the defects, their con- 
centration, wavevector cutoffs, etc. As regards this latter 
quantities, for example, a crude estimate can be obtained 
from the smallest possible wavevector in a system with 
a finite size L: q^-m ~ 27r/L In the case of out-of- 
plane phonons, however, the combined effect of fluctua- 
tions and anharmonicity yields a strong renormalization 
of the bending rigidity k 0, 2^, 2^ . This eventually sup- 
presses out-of-plane fluctuations, making it possible the 
low temperature flat phase of the sheet. This result can 
be seen as due to an anharmonicity-induced infrared cut- 
off , which is the relevant one for large enough samples 
[2^ (i.e., in the thermodynamic limit). But in the case 
of in-planc phonons the tendency is just the contrary. In 
this case, nonlinear effects make the divergence of local 
in-plane fluctuations even stronger In consequence, 
by employing Eq. (fTUl) with Tq hcqmin, one is actually 
underestimating the Koshino- Taylor effect if the size of 
the sample is large enough. 

At this point, it is worth mentioning that a similar de- 
pendence on the temperature is obtained from the renor- 
malization of the impurity scattering due to the Friedel 
oscillations of the corresponding exchange field • But 
in contrast to the Koshino- Taylor effect, such a renormal- 
ization is not operative in the case of Coulomb scatterers 
and does not scale with the size of the system. In any 
case, it is worth noticing there is a correlation between 
the temperature-dependent part of these contributions 
and the residual T = resistivity which might be useful 
in their experimental identification [20| . 

In summary, we have discussed the peculiarities of the 
Koshino- Taylor effect in graphene. Both in-plane dis- 
placements of the defects and deformations of the scat- 
tering potentials due to out-of-plane distortions of the 
graphene sheet yield a T In T dependence of the defect- 
induced resistivity. These two contributions, however, 
differ in their dependence on external parameters such as 
the carrier concentration and the tension of the graphene 
sheet. The key ingredient in the Koshino- Taylor effect is 
the local fluctuations of the lattice which, in graphene, 
can be quite large due to the two-dimensional character 
of the system. Graphene thus represents a model system 



to study the Koshino- Taylor effect, an effect somewhat 
elusive in ordinary metals. 
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